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ABSTRACT 

The transfer of energy and angular momentum in the magnetic coupling (MC) of a 
rotating black hole (BH) with its surrounding accretion disc is discussed based on a 
mapping relation derived by considering the conservation of magnetic flux with two 
basic assumptions: (i)the magnetic field on the horizon is constant, (ii) the magnetic 
field on the disc surface varies as a power law with the radial coordinate of the disc. 
The following results are obtained: (i) the transfer direction of energy and angular 
momentum between the BH and the disc depends on the position of a co-rotation 
radius relative to the MC region on the disc, which is eventually determined by the 
BH spin; (ii) the evolution characteristics of a rotating BH in the MC process without 
disc accretion are depicted in a parameter space, and a series of values of the BH spin 
are given to indicate the evolution characteristics; (iii) the efficiency of converting 
accreted mass into radiation energy of a BH-disc system is discussed by considering 
coexistence of disc accretion and the MC process; (iv) the MC effects on disc radiation 
and emissivity index are discussed and it is concluded that they are consistent with the 
recent XMM-Newton observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15 
with reference to a variety of parameters of the BH-disc system. 
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1 INTRODUCTION 

Recently, the magnetic coupling of a rotating black hole 
(BH) with its surrounding disc has been investigated by cer- 
tain authors (Blandford 1999; Li 2000, Li 2002a, b, hereafter 
Li02a and Li02b, respectively; Li & Paczynski 2000; Wang, 
Xiao & Lei 2002, hereafter WXL). This coupling can be 
regarded a variation of the Blandford-Znajek (BZ) process, 
proposed two decades ago (Blandford & Znajek 1977) . With 
closed magnetic field lines connecting a rotating BH with the 
disc, energy and angular momentum can be transferred from 
the BH to the disc and henceforth this energy mechanism 
is referred to as the magnetic coupling (MC) process. The 
load in the MC process is the surrounding disc, which is 
much better understood than the remote astrophysical load 
in the BZ process. Energy and angular momentum are al- 
ways transferred from the BH to the unknown remote load 
in the BZ process, while the transfer direction of energy and 
angular momentum in the MC process depends not only on 
the angular velocity of the BH but also on that of the disc 
where each closed field line penetrates. 

Usually, the transfer of energy and angular momentum 
in the MC process is stated as follows: 'If a BH rotates faster 
than its surrounding disc, it exerts a torque on the disc and 



energy and angular momentum are extracted from the BH 
and transferred to the disc, and vice versa'. However, this 
statement is rather vague and only applicable to the case 
that the closed magnetic field lines are attached to the inner 
edge of the disc. In a more realistic model the closed field 
lines are assumed to connect the BH horizon with the disc by 
attaching a ring with inner and outer boundary as shown in 
Fig.l, where ri„ and r out are the radii of the inner and outer 
boundary of the MC region, respectively, and 6\ and 82 are 
the corresponding angular coordinates on the BH horizon. 

It is shown in WXL that the width of the MC region 
should not be given at random, which can be determined by 
conservation of magnetic flux in the context of general rel- 
ativity. In this paper the mapping relation derived in WXL 
between the angular coordinate on the BH horizon and the 
radial coordinate on the disc is modified by a boundary con- 
dition on the closed field lines connecting the two loops in 
the equatorial plane of the BH. It turns out that the trans- 
fer direction of energy and angular momentum is determined 
by the position of the co-rotation radius relative to the MC 
region, and ultimately depends on the BH spin eventually. 
In addition, other MC effects, such as BH entropy change, 
efficiency of releasing energy, disc radiation and emissivity, 
also appear to be related intimately to the BH spin . In 
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Figure 1. The poloidal magnetic field connecting a rotating BH 
with its surrounding disc 



particular, as argued in Li02b, the MC effects on emissivity 
are consistent with the recent XMM-Newton observation of 
the nearby bright Seyfert 1 galaxy MCG-6-30-15. We find 
that a variety of MC parameters are suitable for the above 
observation. 

This paper is organized as follows. In Section 2 we give 
the basic evolution equation of a BH in MC process, and 
derive a modified mapping relation between the angular co- 
ordinate on the BH horizon and the radial coordinate on the 
disc. In Section 3 the evolution characteristics of a rotating 
BH in the MC process without disc accretion are depicted 
in terms of the BH spin in a parameter space, and a series 
of values of the BH spin are given to indicate the evolu- 
tion features. The sequence of these values are related to 
the transfer of energy and angular momentum between the 
BH and the disc, and guaranteed by the second law of BH 
thermodynamics (Wald 1984). In Section 4 the efficiency 
of BH-disc systems in converting accreted mass into radia- 
tion energy is discussed in terms of the coexistence of disc 
accretion and the MC process. In addition, we discuss MC 
effects on disc radiation and emissivity index that are consis- 
tent with the recent XMM-Newton observation of the nearby 
bright Seyfert 1 galaxy MCG-6-30-15. Finally in Section 5, 
we summarize our main results. 

In order to facilitate the discussion of the correlation of 
BH spin with MC effects in an analytic way, we make the 
following assumptions: 

(i) The disc is both stable and perfectly conducting, and 
the closed magnetic field lines are frozen in the disc. 

(ii) The disc is thin and Keplerian, and lies in the equa- 
torial plane of the BH with the inner boundary being at the 
marginally stable orbit. 

(iii) The magnetic field is assumed to be weak, and the 
effect of instabilities of the disc and the magnetic field are 
ignored; The magnetosphere is stationary, axisymmetric and 
force-free outside the BH and the disc;. 

(iv) The magnetic field is assumed to be constant on 
the horizon, and to vary as a power law with the radial 
coordinate of the disc. 

Throughout this paper the geometric units G = c = 1 
are used. 



2 BASIC EQUATIONS FOR BH EVOLUTION 
AND MAPPING RELATION IN THE MC 
PROCESS 

As argued by Macdonald & Thorne (1982, hereafter MT), 
the magnetic field on the horizon is brought and held by the 
surrounding magnetized disc. So both disc accretion and the 
MC process should be taken into account in our model for 
BH evolution. Based on conservation of energy and angular 
momentum, the basic equations for BH evolution with the 
coexistence of disc accretion and the MC process are written 
as 



dM/dt = E ms Mn - Pmc 
dJ/dt = LmsMo — Tmc- 



(1) 
(2) 



Combining equations and @ we have the evolution of 
BH spin: 



da,/dt = M- 2 (L ms M D - Tmc) 

—2M~ 1 a,(E ms Mn - Pmc). 



(3) 



Now we give a brief explanation for the quantities in 
equations Q and M and J are mass and angular mo- 
mentum of the BH, respectively, a* = J / M 2 is the dimen- 
sionless angular momentum of the BH, and referred to as 
the BH spin. E ma and L ms are specific energy and angu- 
lar momentum corresponding to the inner edge of the thin 
disc, i.e. the last stable circular orbit, respectively (Novikov 
& Thorne 1973). Md is the accretion rate. At first sight, 
Mo will be affected by the MC process owing to angular 
momentum transferred between the rotating BH and the 
disc. We derived the MC correction on Md by considering 
the conservation law of angular momentum and the angular 
momentum transferred at the inner edge of the disc (see Ap- 
pendix A in WXL) . Unfortunately, the modification to Md , 
equation 142> in WXL, is incorrect, since accretion rate can- 
not depend on radius in a stationary flow, and the validity of 
the following results of sections 5 and 6 in WXL is doubtful. 
In this paper the MC effects on the accretion rate are not 
considered for the following reasons. 

(i) A constant accretion rate is required by the conser- 
vation of mass in a stationary disc. 

(ii) From equations (1611 and (I62H we find that the vis- 
cous torque related to disc radiation can regulates itself ev- 
erywhere to counteract the effects of angular momentum 
transferred into the disc on the accretion rate. 

In the basic evolution equations (0 and @, Pmc and 
Tmc are the rates of extracting energy and angular momen- 
tum from the rotating BH by the MC process, and hence- 
forth are referred to as the MC power and the MC torque, 
respectively. We have derived the expressions for Pmc and 
Tmc in WXL as follows: 



Pmc/Po = 2at 



ir/2 



/3(l-/9) sin 3 Odd 



i 2-(l-<2)sin 2 



Tmc /T = Aa* {I + q) 



r/2 



9l 



(1-/3) sin 3 Odd 
2 - (1 - q) sin 2 6 



(4) 



(•>) 



where 6 is the angular coordinate on the horizon varying 
from 6\ to 7r/2. Pq and To are defined as 
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Pa = (B%) M 2 « B\Ml x 6.59 x 10 44 erg • s"\ 
To = (B%) M 3 w B|M| x 3.26 x 10 47 g ■ cm 2 ■ s" 



(6) 



and B4 and Ms are y and M in the units of 10 4 gauss 

and 1O 8 M0, respectively. The parameter /3 = Qd/£Ih is 
denned as the ratio of the angular velocity on the thin disc to 
that on the horizon, and we have the following expressions: 



n H = a f /(2r H ), r H =M(l+q), q=y/l-al, (7) 
and 

Q.F = = .,, ,, r, 



(8) 



where r H is the radius of the horizon of a Kerr BH. Qf is 
the angular velocity of the closed field line connecting the 
BH and the disc, and flp — firj arises from the freezing- in 
condition in the disc (MT). 

In order to calculate Pmc and Tmc we should first de- 
termine the mapping relation between the BH horizon and 
the disc. Considering the flux tube consisting of two adja- 
cent magnetic surfaces "i" and "i + 1" as shown in Fig.l, we 
have A^>h = A^d required by continuum of magnetic flux, 
i.e. 

B ± 2n {zup) r=rH d9 = -B z 2n (zup/VA) e=n/2 dr, (9) 

where B± and B z are the normal components of magnetic 
field at the horizon and the disc, respectively, and 



(wp) r = (Est 



= 2Mr„ sin0, 



( OT p/VA) =E/VA = a- 1 p, 



where 



a = (pVA/E 



9=ir/2 



1 - 2 Xm lt- 1 + afxgjilZ 
1 + alxTnU- 2 + 2a 2 x^st 



(10) 
(11) 

(12) 



In equation 1121 £ = r/r ms is a dimensionless radial param- 
eter defined in terms of the radius r ms of the last stable 
circular orbit. Following Blandford (1976) we assume that 
B z varies as 



(13) 



Owing to the lack of knowledge of configuration of the mag- 
netic field in the gap region between the horizon and the 
inner edge of the disc, we proposed the following assump- 
tions in WXL: (a) the inner boundary of the MC region is 
located at the inner edge of the disc, with which the closed 
field lines connect the horizon at 9 — n/2; (b) the strength 
of the magnetic field at the horizon is equal to that at the in- 
ner edge of the disc. However the former is probably greater 
than the latter by numerical simulation in the disc (Ghosh 
& Abramowicz 1997, and the references therein). Consider- 
ing conservation of magnetic flux in the inner boundary of 
the MC region, we propose the following relation to replace 
assumption (b): 



2tyt h B± = 2-kvjo (r m<1 ) B z (r„ 



(14) 



where vjd {r ms ) is the cylindrical radius at the inner edge 
of a thin disc and reads 



TUB 



(rms) = MxL y/i + x^ta'l + ^X^a 2 . 



(15) 
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Figure 2. The curves of £ c (dot-dashed line) and £ 

ut versus ci* 

for < a* < 1 with n = 1.1, 1.5 and 3.0 in solid, dashed and 
dotted lines, respectively. 



Equation 1141 implies conservation of magnetic flux corre- 
sponding to the two loops of the same infinitesimal width, 
and the ratio of B± to B z (r ms ) varies from 1.8 to 3 for 
< a* < 1. Incorporating equations 1131 and 1141 we have 

B z =Bx[r H /w D (r in )}C n - (16) 
Combining equations 1101 — 1161 into equation © we have 



sin 9d6 — — G (a, ; £, n) d£, 
where 
G(a*;f,n) 



(17) 



(18) 



2^(l+ai x - 4 a +2ai x - e s )(l-2 x - 



Integrating equation 1171 and setting £ = = 1 at 9 
7r/2, we express the mapping relation by 



<i>m/ = / G(a*;£,n)d£, 



(19) 



and the outer boundary £ ut = r ut/r m s of the MC region 
can be determined by the following equation: 

cos0i=^ G(a»;£,n)d£. (20) 

From equation 1201 we find that £ oui (a*; n, 6%) behaves as a 
non-monotonic function of a* for the given n and 9\ as shown 
in Fig. 2. It is noted that 9\ is taken as 7r/6 throughout this 
paper except where indicated otherwise. 

The co-rotation radius r c is defined as the radius on the 
disc where the angular velocity Hd of the disc is equal to 
the BH angular velocity fijj , i.e. 

2(1 + 9) f 



+ a* 



(21) 



From equation 1211 we can express r c in terms of a 
parameter £ c = r c /r ms as follows: 



£c (a,) = r c /r„ 



-2 
Xms 



(i-«r i/8 (i+«). 



(22) 
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Table 1. Some quantities related to a Rq 



a 


eg 
<I, J 


(out \<C q ) 




R S 


R P 


Hrp 


1.1 


0.2924 


1.28273 


1.10957 


0.63274 


-1.26867 


-1.0104 


1.5 


0.2908 


1.29863 


1.11273 


0.60636 


-1.28244 


-1.0107 


3.0 


0.2835 


1.38780 


1.12758 


0.49025 


-1.35549 


-1.0123 



The parameter £ c decreases monotonically with a, is 
shown by the dot-dashed line in Fig. 2, where there exist 
two kinds of intersections: one is that the curve £ c (a*) in- 
tersects with each curve of £ ut (a*; 0i,n), which is indicated 
by (a° ut ,£out)', the other is that it intersects with £j„ = 1, 
which is indicated by (0.3594, 1). These intersections im- 
ply that the co-rotation radius r c is located at the outer 
and inner boundaries of the MC region for a, — a° ut and 
aT = 0.3594, respectively. The MC region is divided into 
two parts by £ c : the inner MC region (henceforth IMCR) for 
1 < £ < ic and the outer MC region (henceforth OMCR) for 
£ c < £ < £out. Therefore energy and angular momentum are 
always transferred by the closed magnetic field lines from 
the BH into the OMCR with Q.d < Qh, while the transfer 
direction reverses for IMCR with Qd > Qh- The correlation 
of the BH spin with the transfer direction is given as follows: 

(i) For 0.3594 < a, < 1 we have £ c within the inner 
edge, and the MC region is simply the OMCR with the 
transfer direction from the BH to the disc; 

(ii) For a% ut < a, < 0.3594 we have 1 < £ c < £ out , 
and the transfer direction is from the BH to the disc in the 
OMCR, while it is from disc in the IMCR to the BH; 

(iii) For < a* < a° ut we have £ c > £ otlt , and the MC 
region is simply the IMCR with the transfer direction from 
the disc to the BH. 



3 BH EVOLUTION IN THE MC PROCESS 
WITHOUT DISC ACCRETION 

3.1 Equilibrium spin of a BH in MC process 
without disc accretion 

In order to highlight the MC effects on BH evolution we 
discuss a specific case where disc accretion is absent. In this 
case a rotating BH can attain its equilibrium spin a eq in the 
MC process for 1 < £ c < £ out . The basic evolution equations 
and become 



dM/dt = -P MC , 
dJ/dt = -Tmc, 

da,/dt = -M~ 2 T M c + 2M~ 1 a,P M c- 
Setting da*/dt = 0, we have 
Tmc = 2Ma*P M c- 



(23) 
(24) 
(25) 

(26) 



Substituting equations 0, JKJ and 11911 into equation 126H . 
we can obtain at q for the given parameters 9\ and n by 
resolving the following equation: 



s: 



^ (l-P)(l-P + 0g) 
2 esc 2 0-(l-q) 



G(a,;£,n)d£ = 0. 



(27) 



In order to discuss the transportation of energy and angular 
momentum at at q , we define the following ratios: 



r> D* n / T)OUt t~j rriin / rriout 

tip — r MC I r MC , ilT — J MC 1 1 MCi 



(28) 



where Pmc an d Ijfc are the rates of transferring energy and 
angular momentum in the IMCR, respectively, and Pmc 
and Tmc are * ne rates of transferring energy and angular 
momentum in OMCR, respectively. By using the mapping 
relation 1191 these quantities can be expressed as follows: 



r>™ ,2 P / /3(l-/3)G(a»;£,n) 
t^MC = &a.*lro I —z 5-7 — t. r~ 



mc = 2a, P 



2 esc 2 — (1 - q) 

/3(l-0)G(a.;Z,n) 
2csc 2 6>- (1-g) 



Tmc = 4a*T (1 + q) 



(l-/7)G(a.;g,n) 
2 esc 2 6 — (1 - q) 



(29) 
(30) 
(31) 



T MC - 4a,T (1 + q) 2csc 2 e _ (1 _ g) d t- (32) 

From equations G3, ll27t— l|3ljl we obtain some quan- 
tities related to af as listed in Table 1. 

From Table 1 we obtain the following results: 

(i) Both Rp and Rt are less than —1, and these imply 



Pmc (al q ) = Pll c (a?) + Pgfg (a?) < 0, 
Tmc (at q ) = T™ c (of) + TJEfg (at q ) < 0, 



(33) 



which means that the transportation of energy and angular 
momentum in OMCR is dominated by that in IMCR, i.e. 
they are transferred as a whole from the disc to the BH at 
equilibrium spin a% q . 

(ii) Defining R^ = (f c — l)/(£ «t — Cc) to indicate the 
ratio of the radial width of IMCR to that of OMCR, we find 
that both R^ and a% q decrease as the increasing n, while 
both £out (a* q ) and £ c (a% q ) increase with increasing n. These 
results are related to the conservation of magnetic flux with 
the fact that more magnetic field is concentrated in IMCR 
for greater value of n. 



3.2 BH evolution in the MC process without disc 
accretion 

We can discuss the BH evolution in the corresponding pa- 
rameter space as proposed in WXL. Evolution equations 
1231 1. 1- It and 1251 can be rewritten as 



dM/dt = -Pmc = 2P f (a,; n, Oi) , 



(34) 
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11 

Figure 3. Parameter space for BH evolution in the MC process 
with 

1.1 <n < 2.1 and 0.2881 < a* < 0.2960. 

Table 2. The signs ofdM/dt, dj/dt, and da,/dt 
in the regions of parameter space 



in the IMCR or in the OMCR. In fact, as argued above, the 
transfer direction of energy is always the same as that of 
angular momentum either in the IMCR or in the OMCR, 
and the opposite signs of Pmc and Tmc rest in the fact 
that Pfifc and T^ic are dominated by P M l c and T£g for 
the above value range of the BH spin, respectively. We can 
explain this order of the BH spin by using the second law of 
BH thermodynamics in the next section. The order of the 
above specific values of BH spin is given as follows: 



r, , out - en „ T „ P - i 

< a, < a, < a t < a t < a, 



< 1. 



(37) 



We can use the parameter space to determine the evo- 
lution characteristics of a rotating BH in the MC process, 
provided that the initial BH spin a* and the power law index 
n are given. For example a BH with initial spin af < a» < 1 
will evolve to a% q and reach the curve g = eventually by 
passing region IV, III and II one after another. So we can de- 
termine its evolution characteristics in each evolution stage 
represented by the corresponding region in Fig. 3. 

The status of MC region and the signs of Pmc and Tmc 
for different value range of a* are shown in Table 3. 



Region 



Position 



da,/dt dj/dt dM/dt 



I left of g = > > > 

II between g = and h = < >0 >0 

III between h = and / = 0<0 <0 >0 

IV right of/ = <0 <0 <0 



dJ/dt = -Tmc = 4T h (a* ; n, 0i) , 

h (a,;n, 01 ) = -a, (1 + q) ^3f-^ ^, 

da*/dt = AT M~ 2 g (a,;n,6±) , 



g (a*;n,0i) = -a„ (1 + q) X 

[■tout (l-f3)(l-P+l3q)G(a,;a, n ) 
Jl 2 csc^ S-(l-q) 



(35) 



(36) 



From equations 1341 — 136H . we find that /(a*;n, 9i), 
h(a,;n,9i) and g(a*\n,6\) have the same signs as 
those of dM/dt, dj/dt and da*/dt, respectively. Setting 
/ (a»;n,#i) = 0, h(a,;n,6i) — and <? (a*; n, #i) = for 
#i = 7r/6, we have three characteristic curves in the two- 
dimension space consisting of the parameters a» and n as 
shown in Fig. 3, where each black dot with an arrowhead 
represents one BH evolution state. From left to right, the 
parameter space is divided into four regions by these three 
curves, and the signs of the rates of change of M, J and a, 
in these four regions are listed in Table 2. 



From Fig. 3 we find that Pmc changes its sign from 
negative to positive at af (/ = 0) and Tmc does it at 
(h — 0). The inequality, < af , holds for all possi- 
ble values of n, because the curve h — is located on the 
left of / = as shown in Fig. 3. This result means that 



Pmc is negative with the positive Tmc in the value range 
< a» < . However, it does not mean that the transfer 
direction of energy is opposite to that of angular momentum 



3.3 BH entropy change in the MC process 
without disc accretion 

Entropy Sh of a Kerr BH can be expressed as (Wald 1984; 
Thorne, Price & Macdonald 1986) 



Sh = 2ttM 2 (1 + q). 



(38) 



From equations 1231 and 12411 we derive the rate of change 
of Sh as follows: 

dSn/dt = e H x (QhTmc - Pmc) 
= 87rMP (l+g) 2 (q- 1 



8 2 (1-/3) 2 sin 3 Bd6 
2-(l-q)sin a ' 



(39) 



where Oh 



1)C 

is the temperature on the BH hori- 



47rA<f(l + (j) 

zon. Inspecting equation l|39^ we find that dSn/dt > 
always holds. The contribution to dSn/dt arises in two 
parts: one part from Tmc, which is positive and negative 
for a% < a, < 1 and < a» < aif, respectively, and the 
other part from —Pmc, which is positive and negative for 
< a» < and a* < a* < 1. It is noticed that the two 
contributions are all positive for < a* < . So the or- 
der a% < a* < af is guaranteed by the second law of BH 
thermodynamics, otherwise the reverse order will result in 
dSn/dt < 0, and the second law of BH thermodynamics will 
be violated. The dimensionless rate of change of Sh can be 
written as 



(rate) , 



dS H /dt 
(dS„/dt) 



(i + q) (q - 1 )Ji 2c E cifl-(i-,) d ?> 



(40) 



where (dS H /dt) = SttMPq » 1.07 x 10 60 Bf Af 8 3 er fl ■ K' 1 ■ 
s _1 . From equation 14U1 we obtain the curves of (rate) s 
versus a* for different values of n as shown in Fig. 4. 

From Fig. 4 we find that (rate) s varies non- 
monotonically as a*, attaining its minimum (rate)™ 111 at 
a™ m and maximum (rate)™ ax at a™ 1 " 1 , respectively. These 
limit values corresponding to the different values of n are 
listed in Table 4. 
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Table 3. The signs of Pmc > ^MC an d the rates of change of BH parameters. 
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Figure 4. The curves of (rate) s versus a, with n = 1.1, 1.5 and 
3.0 in solid, dashed and dotted lines, respectively, with the upper 
panel for < a* < 1, and the lower panel for 0.99 < a» < 1. 



Table 4. The values of (rate)™™ , (rate) max 
and the corresponding BH spin 



n 




(rate) min 


{p 1 * )max 


(rate) max 


1.1 


0.2909 


0.0005382 


0.99959 


0.44995 


1.5 


0.2892 


0.0005794 


0.99960 


0.47772 


3.0 


0.2812 


0.0008010 


0.99887 


0.71611 



From Table 4 we find that (af ) is located between 

V / min 

a° ut and af , which corresponds to increasing M , J and a* 
with Pmc < and Tmc < 0, while (af ) is very close 

v / max 

to unity, corresponding to the decreasing M, J and a* with 
Pmc > and Tmc > 0. From Fig. 4 we notice that (rate) s 
decreases rapidly when a„ > (a?) , and we shall explain 
this result in Section 5. 



4 MC EFFECTS ON THE BLACK HOLE-DISC 
SYSTEM 

Now we are going to discuss some evolution characteristics 
in MC process with disc accretion (henceforth MCDA). The 
interaction between a rotating BH with a disc is crucial for 
at least two reasons: 

(i) The magnetic field on the BH horizon is brought out 
and maintained by the surrounding magnetized disc; 

(ii) The transfer of energy and angular momentum be- 
tween the BH and the disc will remarkably affect not only 
BH evolution but also disc radiation. 

As the magnetic field on the BH is supported by the 
surrounding disc, there are some relations between Bh and 
Md- As a matter of fact these relations might be rather 
complicated, and would be very different in different situ- 
ations. One of them is given to investigate the correlation 
between BH spin and dichotomy of quasars by considering 
the balance between the pressure of the magnetic field on 
the horizon and the ram pressure of the innermost parts of 
an accretion flow (Moderski, Sikora & Lasota 1997), i.e., 

B 2 H /(8tt) = Pram ~ pc ~ M D / (4^%) . (41) 

From equation (1411 we assume the relation as 



M D = B 2 H M 2 (l+qf /2. 



(42) 



Substituting equations J1J , J5J and 1 1-1) into equations ijTjl 
we have 



(43) 
(44) 
(45) 

(a*,n,0i) and 
g MCDA (a*, n, 8i) are the characteristics functions of 
BH evolution in the MCDA. 



dM/dt = f MCDA (a*,n,ei)M D , 
dJ/dt = Mh MCDA (a,,n, 0i) M D , 
da*/dt = M~ 1 g MODA (a*, n, 9i) M D , 
where f MCDA (a f ,n,6 L ), h M 
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Figure 5. The curves of r\ DA (solid line) and r\ MC versus a* with n = 1.1, 1.5 and 3.0 in dashed, dot-dashed and dotted lines, 
respectively. (a) < a, < 1, (b) 0.288 < a» < 0.297, (c) 0.895 < a» < 0.900, (d) 0.988 < a» < 0.999 



Table 5-Some specific value ranges of r\ MC and corresponding values of a* 



Value range of a* 



Vmc < 



< V M c < Vdj 



(.VMc) max { al l)r, 



1.1 (0,0.2960) (0.2960,0.8977), (0.9966,1) (0.8977,0.9966) 0.3070 0.9949 

1.5 (0,0.2947) (0.2947,0.8974), (0.9969,1) (0.8974,0.9969) 0.3101 0.9951 

3.0 (0,0.2890) (0.2890,0.9002), (0.9976,1) (0.9002,0.9976) 0.3182 0.9962 



4.1 Efficiency of the BH-disc system 

The total efficiency of converting accreted mass into radi- 
ation energy corresponding to the BH evolution in MCDA 
is 



ris = 1- (dM/dt)/M D = 1 - f MCDA (a*, Mi) , 



(46) 



and it consists of two parts due to the MC process and disc 
accretion as follows: 



1)* = Vda + Vmc 
where 

Vda = 1 - E ms, Vmc = Pmc/ M l 



(47) 



(48) 



Substituting equations Q an d tl 1 — D into equation l|48ll . we 



express rj A 



(a*,n) = 



4(1-9) r ut 0(l-/3)G(<u;t,n) 



(1 + 9) 



2csc 2 6»-(l-g) 



df. (49) 



ft is noticed that rj MC only depends on the parameters o* 
and n in our model. The curves of r/ DA and r\ MC versus a* 
for the given n and 9\ are shown in Fig. 5. 

From Fig. 5 we find the following results for the efficiency 
in the MCDA: 

(i) The efficiency r\ DA increases monotonically with a,, 
while rj MC varies non- monotonically with a,, attaining a 
maximum for a high BH spin. 

(ii) The efficiency rj MC , is greater than r\ DA for some 
value range of the high BH spin as listed in Table 5. 

(iii) r? MO = occurs just at af for Pmc = 0, and rj MC < 
occurs for < a» < af , implying that energy is transferred 
from the disc to the BH by the magnetic field. 

(iV)The efficiency r\ MC decreases very rapidly from its 
maximum to zero after a* passes across the maximum point. 



(v) Although rj uc is not sensitive to the variation of the 
parameter n, the values of (t| MC ) m „ and the corresponding 
(a2) become greater as n increases. 

V * / max o 

Some specific value ranges of r\ MC and the correspond- 
ing a, are listed in Table 5. 



4.2 MC effects on disc radiation and transfer of 
angular momentum 

Based on the three conservation laws of mass, energy and 
angular momentum the following equation of radiation from 
a thin disc was derived by considering the MC effects in 
Li02a: 



F = F DA + Fmc, 



(50) 



where Fda is the radiation flux due to disc accretion as 
given by Page and Thorne (1974, hereafter PT74): 



F DA = M D f DA /(4 7IT ), 



(51) 



and Fmc is the radiation flux due to the MC effects and 
expressed by 

, -2 



X 



Fmc = (E+ -fl D L 

fl (E+ -n D L+) Hrdr 



(52) 



Function H in equation i1o2l is the flux of angular momen- 
tum transferred between the BH and the disc by the mag- 
netic field. E + and L + are the specific energy and angular 
momentum of a particle in the disc, respectively, and read 
( Novikov & Thorne 1973) 



E^ 



l-2 X - 2 + a« X " 3 
(l-3x" 2 +2a,x _3 ) 1/2 ' 



(53) 
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T + M x (i-2a,x~ 3 + »*x" 4 ) 

(l-3x- 2 + 2a»x- 3 ) 1/2 



(54) 



and we have E + — E ma and L + = L m3 for f = 1 with 

X = X">»- 

Very recently Li pointed out that the magnetic cou- 
pling between a black hole and a disc can produce a very 
steep emissivity with index a = 4.3 ~ 5.0, which is consis- 
tent with the recent XMM-Newton observation of the nearby 
bright Seyfert 1 galaxy MCG-6-30-15 (Li02b). However such 
a steep emissivity is very difficult to be explained by a stan- 
dard accretion. The emissivity index is defined as 



a = —dlnF/dkir, 



(55) 



which mimics F oc r~ a locally. In Li02b the calculation for 
the emissivity index is done for a stable non-accretion disc, 
and the flux function H is assumed to be distributed from 
v = r ms to r = r out with a power law: 



H = 



Ar n 



0, 



< r < r oll 



(56) 



r > Ton 



where A is regarded as a constant in Li02b. 

Equation I56H can be modified by using the mapping 
relation in our model. From the conservation of angular mo- 
mentum and equation JSJ we have 

^/fr^^-MMM^?, (57) 



where 

36/ dr = (39/30 (d£/dr) 



2 - (1 - q) sin^ 
G(o.;£,n) 



dr 



(58) 



Substituting equations 11811 and 15811 into equation 1571 . we 
have 



H(a*;Z,n)/H 







1 < £ < £o«i 



(59) 



where we have H = (B 2 H ) M = 1.48 X 10 21 x B\M % g-s~ 2 , 
and 



a,(l-/3)(l + 9 ) 
2irx^ s [2csc 2 e-(l-ij)J 



(60) 



8,-3 



V / (l + ^X^+2^X^)(l-2x^ 2 s «- 1 +^X^|- 2 )' 



Since the magnetic field decreases with increasing r, the 
power-law index in equation 1561 should be negative, obey- 
ing m = —n. Thus we derive a modified expression 1591 for 
the function //, where the coefficient A(a*,£) is dependent 
on a* and £ rather than a constant. For the given values of 
o* and n we have the curves of H (a,; £, n)/Ho versus the 
radial coordinate £ for 1 < £, < £ ou t as shown in Fig. 6. 

Inspecting Fig. 6, we have the following results on the 
transfer of angular momentum: 

(i) The value of H (a* ;£,n)/Ho is negative for very low 
BH spin as shown in Fig. 6(a), and its value is positive for a 
very high BH spin as shown in Fig. 6(c). 

(ii) As shown in Fig. 6(b), positive H (a„; £, n)/Ho co- 
exists with negative H(a*;£,n)/Ho for BH spin, such as 
a» = 0.328, which is consistent with the value for the co- 
existence of the OMCR and the IMCR as shown in Table 
3. 

(iii) As shown in Fig. 6b and 6c, the flux of angular 
momentum transferred decreases as £ approaches £ ut- 

By using equations 151L 1521 and 159H we have the 
curves of lg(F M c/F ) and lg(F DA /F ) versus lg(r/M) as 
shown in Fig. 7, where Fo is defined as Fo = Mn/r^s. 

Combining equations 152L 1551 and I59L we have the 
curves of the emissivity index a versus lg (r/M) as shown 
in Fig.8. 

From Fig.7 we find that lg (Fmc/Fq) varies much more 
steeply with lg (r/M) than lg (Fda/Fq) does. It is found in 
Fig.8 that we can adapt the emissivity index arising from 
— d In Futc/d In r to the observations of the observations by 
the curves in the shaded region, while the index of a standard 
accretion disc (SAD) is far below the shaded region. This re- 
sult is one of the observational signatures of the magnetic 
coupling between a rotating BH and its surrounding disc. 
Compared with the model given in Li02b we have more pa- 
rameters to choose for adapting the emissivity index to the 
observations. Considering both the mapping relation 1191 
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Table 6. The MC parameters adapting the emissivity index to the observations 



81 


a* = 0.990 


a* = 0.992 


a» = 0.994 


a* = 0.996 


a* = 0.998 


jr/12 


2.340 


2.329 


2.326 


2.343 


2.423 


tt/4 


3.335 


3.339 


3.359 


3.413 


3.576 


2tt/5 


7.563 


7.673 


7.851 


8.176 


8.944 



a , = . 995 



V. 

SAD \\\ 




3.2 0.3 0.4 

lg{rlM] 



(a) 




lg{rlM\ 
00 

Figure 7. The curves of 1%(Fmc/Fq) versus lg(r/M) with £ ou t = 
2.5 and n = 6.5, 7 and 7.5 in dotted , dashed and dotted-dashed 
lines, respectively. The curve of lg(Frjj^/ Fo) versus lg(r/M) is 
plotted by a solid line indicating a standard accretion disc (SAD): 
(upper panel) a* = 0.995, (lower panel) a* = 0.998. 



and the observations of the Seyfert 1 galaxy MCG-6-30-15, 
we have the values of the power-law index n corresponding 
to the different values of 81 and a* as shown in Table 6. 

From Table 6 we find that the values of n adapted to 
the observations are generally less than those given in Li02b, 
if 0i is not close to n/2. This result implies that the concen- 
tration of the magnetic field on the disc could be relaxed in 
our model. 

It has been proved in PT74 that the internal viscous 
torque per unit circumference is related to the radiation 
flux F by 



2 (E~ l 



n D L+) 



(-dQo/dr) 



(61) 







fl.=0 . 995 


MCG-6-30-15 

s 




//.■■ 

■ /■■ 
It: 
■/■■ 




II: 
■/■' 
//;' 

f 

i 




- f 
t 




0.2 


3 °' 4 0.5 
lg[ rIM ] 

(a) 




Figure 8. The emissivity index a versus lg(r/M) with t; ut = 
2.5 and n = 6, 7 and 8 in dotted, dashed and dot-dashed lines, 
respectively. The emissivity index of the Seyfert 1 galaxy MCG-6- 
30-15 inferred from the observation of XMM-Newton is shown by 
the shaded region. The emissivity index of a standard accretion 
disc (SAD) is plotted in solid line, (upper panel) a* = 0.995 , 
(lower panel) a„ = 0.998. 



Equation 1611 was derived based on the three laws of con- 
servation, and it is proved to be valid for the MC process in 
Li02a. Combining equations I5U1 and 1611 . we can express 



the contribution to WZ as follows: 



w; = {w;) 1 



MC ' 



(62) 



where ( WZ) and (W^) are the contribution due to 
disc accretion and the MC process, and are related to Fda 
and Fmc by equation 1611 . respectively. Since the flux of 
angular momentum transferred away from the disc by radi- 
ation (henceforth FAMFD) is FmcL + , and that transferred 
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5 « 6 
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Figure 9. The curves of Rp versus £ for 1 < § < £ ut with n = 1.1, 1.5 and 3.0 in solid, dashed and dotted lines, respectively, (a) 
a, = 0.1, (b) a, = 0.328, (c) a, = 0.998. 



from the BH into the disc by the magnetic field (henceforth 
FAMFH) is H, the ratio of the two can be expressed as 



(63) 

= -§tfe ( S + _ n D L+y 2 /; ms (£+ - « D L+) Hrdr. 

From equations 1631 we obtain the curves of Rf versus £ for 
the given values of n and a, as shown in Fig. 9. 

Combining Figs. 6 and 9, we obtain the following rela- 
tions between FAMFD and FAMFH: 

(i) As shown in Fig. 9(a), we have < Rf < 1, which 
means that the absolute value of FAMFD is less than that 
of FAMFH, i.e. H < F M cL + < 0, for a BH with low spin 
such as a* = 0.1. 

(ii) As shown in Fig. 9(b), we have < Rf < +oo and 
— oo < Rf < 0, and Rf approaches infinity at the place 
where H changes its sign. So we infer that FmcL + < 
holds in the whole MC region 1 < £ < £ ou t by considering 
the sign of H in Fig. 6(b). 

(iii) As shown in Fig. 9(c), we have < Rf < 1 for 
& n < £ < £e 9 , where FAMFD is dominated by FAMFH. 
And we have R F > 1 for £ eq < £ < £ out , where FAMFD 
dominates over FAMFH. It is easy to obtain that the radial 
coordinate £ e<3 indicating Rf = 1 is equal to 1.498, 1.552 
and 1.857 for n =1.1, 1.5 and 3.0, respectively. 

From the above discussion we find that both FAMFD 
and FAMFH are related intimately not only to the BH spin 
but also to the disc location. 



5 SUMMARY 

In this paper the transfer of energy and angular momentum 
between a rotating BH and its surrounding accretion disc is 
discussed in detail by considering MC effects. Our discus- 
sion is given for two cases: (i)the MC process without disc 
accretion, (ii)the MC process with disc accretion. In the first 
case only the two conservation laws of energy and angular 
momentum are used, while in the second case the conser- 
vation of accreted mass is used in addition. Compared with 
Li02a and Li02b, the mapping relation (1191 is used to depict 



the MC process throughout this paper, and the correlation 
of some parameters with MC effects is discussed, where the 
BH spin o* , the power-law index n and the radial coordinate 
£ are involved. 

Compared with the BZ process we find that the MC 
effects do not increase monotonically as the BH spin. For 
example, in Figs. 4 and 5 both (rate) s and r\ MC attain their 
maxima and then decrease very rapidly as a* approaches 
unity. These results can be explained by the equations 1401 
and 1491 with the behavior of £ otl t as a* approaches unity. 
From Fig. 2 we find that the outer radius of MC region 
approaches the inner edge of the disc very closely as a* 
approaches unity, and accordingly the ratio /3 = Qd / 
approaches unity along as r ms approaches the horizon 
radius r H . 
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